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We analyze the nonequilibrium transport properties of a parallel double quantum dot in terms 
of its full counting statistics (FCS). The parameters of the setup are assumed to be such that both 
subsystems are driven into the Kondo regime. After a series of transformations the Hamiltonian is 
then mapped onto a Majorana resonant level model, which effectively describes the Toulouse point 
of the respective double impurity two-terminal Kondo model. Its FCS is then obtained at arbitrary 
constellation of voltage, temperature, and local magnetic fields. We identify two different transport 
processes corresponding to single electron tunneling as well as an electron pair process and give 
the respective effective transport coefficients. In the most universal linear response regime the FCS 
turns out to be of a binomial shape with an effective transmission coefficient. Furthermore, we 
find a complete transport suppression (antiresonance) at a certain parameter constellation, which is 
similar to the one found in the noninteracting quantum dots. By an explicit expansion around the 
Toulouse point we show that the antiresonance is universal and should be observable in the generic 
Kondo dot setup. We discuss experimental implications of our predictions as well as possible routes 
for generalizations of our approach. 

PACS numbers: 73.23.-b, 73.63.Kv, 85.75.-d 



I. INTRODUCTION 

Quantum dot (QD) is an important paradigm in many 
research fields. One of the most important application 
areas of this concept is the quantum transport. Here the 
QD often represents a basic building element of micro- 
and nanoelectronic circuitrj*^. For the applications it is 
not only important to understand the physics of individ- 
ual devices but also of more complicated arrangements 
of QDs. The focus of our study is the structure with 
the next to single QD complexity level - which is besides 
slight modifications and similar geometries^"— a double 
QD structure in parallel geometry.—"— 

The most widespread modeling strategy for a QD is 
the bottom-up approach. The structure is essentially as- 
sumed to be a zero-dimensional object, which is mod- 
eled by a single spin-degenerate fermionic level. This 
"resonant level" model is then coupled via particle ex- 
change to source and drain electrodes as well as electo- 
statically to a gate electrode which tunes the resonant 
level energy thereby controlling the transmission coeffi- 
cient of the structure. However, because of a confined 
geometry one has to include electron interaction terms 
at least in the minimal form by introduction of finite 
energy cost for the double population of the QD. This 
leads directly to the celebrated Anderson impurity model 
(AIM)^iii which was already applied to parallel double 
dot structures fi^"— It shows up an enormous manifold 
of very different transport regimes, for many of which no 
exact analytical results in nonequilibrium are available. 

Probably one of the most interesting cases is the Kondo 
regime, when the QD is populated by a single electron 
and energetic cost of double population as well as of emp- 
tying the dot are so high that the transport channel in- 
volving only single-electron tunneling events is blockedj^l 



Then only virtual double population is allowed and the 
remaining transport channel is the spin-flip tunneling. It 
has very profound consequences on the transport char- 
acteristics of the systems. For instance, the nonlinear 
I{V) as well as the shot noise turns out to contain two 
different transport channels - single electron tunneling 
and a pair process;^"— which can be seen explicitly in 
their FCSi^i^ The resulting highly nontrivial Fano fac- 
tor 5e/3 has already been verified experimentally Due 
to the additional degrees of freedom of double QDs it 
might, at least in principle, be possible for the electrons 
to be transferred in groups of not only two but also of 
three and even four particles. FCS is most useful in an- 
swering such questions. 

Although the low-energy transport characteristics of 
the Kondo regime of double QDs are by now fairly good 
understood, the nonlinear response and especially the 
FCS of such systems are notoriously difficult to obtain. 
There is, however, a parameter constellation, which is 
on the one hand a non-trivial one in a sense that it cap- 
tures most of the relevant physics of the Kondo limit, and 
which on the other hand allows for a complete analytical 
solution with elementary means. This is known under 
the name of Toulouse point solutions^"—. 

A similar calculation for the double QD setup, be- 
ing an analytic solution, would yield invaluable infor- 
mation about the transport properties of the system as 
well as become an important benchmark for other ap- 
proaches. To the best of our knowledge this kind of cal- 
culation was not yet attempted. There are several ob- 
stacles which need to be circumvented. The first issue 
is the validity of the Hamiltonian in which two differ- 
ent localized spins are simultaneously locally coupled to 
two electrodes. This part of the program was success- 
fully mastered in Ref. [s^. The second question is that 
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FIG. 1: Parallel double quantum dot. The two quantum dots 
are coupled symmetrically (under exchange of the dots) to 
metallic leads. Hi^2 (small arrows) are the in-plane magnetic 
fields. The dot spins are represented by larger arrows. The 
coupling is represented by the doubly headed arrows. 



of the principal applicability of the sequence of trans- 
formations which are usually used in order to obtain an 
essentially quadratic Hamiltonian. As is argued in the 
next section one can indeed apply them after a minor 
adjustment. The last issue is slightly more technical and 
concerns the method for FCS calculation, how one has 
to introduce the counting fields in order to obtain mean- 
ingful results. This is discussed with all necessary de- 
tails in Section ini In Section Hill we present the results 
and discuss the emerging physics. Conclusions Section 
summarizes our findings and gives directions for future 
research. 
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ipl^g{x) is the fermionic field creation operators in the 
electrode a. The bare (without impurity coupling) elec- 
tron dispersion in the leads is linearized around the Fermi 
edge as we are only interested in the low energy sector 
of the system. The interaction term includes the gener- 
alized spin densities s^p{x) = ^il^lsi^)i'^)ss''^l3S'{x)- fJ-B 
is the Bohr's magneton and gimp is the impurity Lande 
factor. We restrict our calculations to the cases where 
the coupling constants fulfill the restrictions 

jap _ TQ/3 _ jaP jLR _ jRL _ r, 

and Jf ^ = ^ = (5) 

first discussed in [s^]. We employ the strategy outlined 
there and proceed by performing a number of transfor- 
mations in order to map the Hamiltonian ([1]) onto the 
one of a noninteracting system. First the Hamiltonian is 
bosonized using the bosonization identity^i^i 



II. MODEL AND CALCULATION METHOD 



The simplest meaningful starting point is the dou- 
ble impurity Anderson model with two electrodes. It is 
shown in Ref. [s^l that by a dedicated Schrieffer-Wolfi' 
transformation^^ the low energy sector is that of two 
magnetic impurities coupled by exchange terms to the 
spin densities in the electrodes (as long as only one elec- 
tron populates the dot). This situation is shown in Fig.[T] 
which also serves to explain the notation conventions. 
We want to emphasize at this point that the applied 
magnetic fields, which are independent on each of the 
dots, are aligned parallel to the plane of the structure 
and therefore do not create an Aharonov-Bohm phase. 
The applied bias voltage V is indicated by the different 
chemical potentials on the leads: ^ f^R — eV. Then 
the corresponding Hamiltonian is given by 



^ — ^kin + ^int + ^mag, (1) 

where the three parts are given by (we use units in which 
e ^ h = vp ^ 1, where vp is the Fermi velocity in the 
electrodes) 




where a is the lattice constant of the underlying lattice 
model. In the next step we introduce the bosonic fields 
4'c,4's,4>f and (corresponding to "charge", "spin", 
"flavour" and "spin-flavour" channels) as linear combi- 
nations of the bosonic fields appearing in 

We apply a unitary transformation which is known to 
change the scaling dimensions of the coupling thus en- 
abling its refermionizatioui^ In the spirit of Ref. [sHi 
we construct it as a product of of two Emery-Kivelson 
rotationsS^ 

[/ = g§x=(^r+^2"). (8) 

around the z-axes of the spins. Here Xs denotes the 
combination of the bosonic field with its corresponding 
Klein factor field Xi/ = 0i^(O) — fiy- One can define new 
fermionic fields 

illidl+4){di+d2) 

iPf = = (9) 

V27ra 

iTT{d{+dl){di+d'2) 

^sf = j= e-'(^=/-^=/) , (10) 
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which allow for a rcfermionization of the Hamiltonian 
yielding 



(-1) 



+ j^^(v4 - + rfi + 4 + f^2) 

+ J"(V'l/ - i^sf){d\ +di+dl + d2) 
+ ^ [AisffiJ, - (J. - 2^) : tPlmMO) ■] {d]d,-l/2), 

3 = 1.2 

(11) 

where the free fermion field ips {x) has to be evaluated at 
X = 0. Now we simplify this expression by introducing 
Majorana fermions for the dot and lead fermions which 
are defined in the following way: 



point by an expansion around the Toulouse point j^i^i^ 
We follow this path in the remainder of the paper. 

At this point we would like to get back to our ear- 
lier restriction of the magnetic fields being aligned to 
the plane of the structure to ensure that no Aharonov- 
Bohm phase (AB phase) is generated. Now, having per- 
formed the bosonization, rotation and rcfermionization 
procedure, we may ask if it is possible to include this 
phase in this formalism. The idea then is to break the 
symmetry J^^ = J^^ and to equip each of the corre- 
sponding terms in the Hamiltonian with an appropriate 
phase e^'"/^. The other terms are not affected by the AB 
phase. Executing the same steps up to the rotation leads 
to an expression that cannot be refermionizcd in the same 
manner as above because of the Klein factors. They lead 
to extra factors of rf 2 in the refermionizcd expression 
and therefore the Hamiltonian cannot be mapped onto a 
nonintcracting system. That is why we restrict ourselves 
to the in-plane magnetic fields. 



For purposes that will become clear later in this work we 
want to define a new parameter 7 to be 7 := — 27r. In 
this new notation the transformed Hamiltonian reads 



H^n - iAiai&i 
iK+risf{bi + ^2) 



- iA2a2b2 + iKLR^fiai 



a2} 



+ 7: ^1(0)^.(0) :(rf-fT|). 



(12) 



We note that after the transformations the coupling con- 
stants Kj as well as the fermionic fields ij^i, have the 
physical dimension ^energy. This can be seen in the 
definitions 



(13) 



where a is the constant that corresponds to the lattice 
spacing that already appeared in Eq. ([6]). 

In its new form the Hamiltonian, which does not 
contain any approximations or simplifications yet, 
possesses only one term which is not quadratic in the 
fermionic fields. In fact, in the case 7 = or J2 = 27r 
it becomes purely quadratic and thus exactly solvable 
by elementary means. This particular point in the 
parameter space is referred to as the Toulouse point— 
Despite its relative simplicity, Toulouse point solution 
carries all features of the generic Kondo effect because 
Jz = 2tt corresponds to rather strong correlations. The 
conventional strategy is to solve the 7 = case first and 
then analyze the robustness of the solution beyond that 



III. FCS OF THE DOUBLE QD SYSTEM 

One of the most generic transport properties of QDs is 
the FCS. It is usually obtained in form of the cumulant 
generating function (CGF) Inx(A), which, being derived 
n times with respect to the counting field A yields the 
nth cumulant (irreducible moment) of the probability to 
measure the transmission of charge Q during a very long 
measurement time 



HQ''' 



1 a" 



■Inx(A). 



(14) 



The current through the constriction is then found from 
^(^) = the shot noise is related to the second 

cumulant at zero temperature S{V) = {{Q^))/T etc. It is 
also useful to construct a Fano factor, which is given by 
the (in our convention dimensionless) ratio of two lowest 
order cumulantsii^ 



F{V) 



SiV) 



(15) 



2IiV) ■ 

It is related to the eminent Schottky formula, which 
allows (at least in principle) a measurement of the 
charge of current carrying excitations, see e. g. Ref. [27l |. 
The main advantage of the CGF is therefore that all 
cumulants are accessible via simple derivation with 
respect to A and setting A = afterwards. 

There are several methods of CGF calculation. We 
shall follow the one presented in Ref. [s^. First the 
counting field A is introduced on the Keldysh contour 



X{t) ^ e{t)9{T ~ t) 



A = A_ , t on forward path 
— A = A+ , t on backward path 

(16) 
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The charge counting operator 



(17) 



is constructed from '7i?(L) parts of the Hamihonian, 
which are responsible for particle transport from the right 
electrode to the left one and vice versa. In the present 
situation it amounts to the decoration of tunneling terms 
V'lcrV'flcr' (V'flcr''/'i<T') in the starting Hamiltonian ([T]) by 
factors 6"*"^^*)/^ (e''*'^*^/^), respectively. 

The CGF can then be derived using the formula^i^ 



x(A) = (Tce-/c^^(*)d*) 



(18) 



where Tc is the Keldysh contour ordering operator and 
the average is taken with respect to the full T-L. Al- 
though this average can be calculated directly, the adi- 
abatic potential method results in a much more com- 
pact algebrai^ It turns out that due to the Feynman- 
Hellmann theorem^ in the limit 7~ — ^ oo there is an 
identity Inx(A) = —iTU{\-,\+), where the adiabatic 
potential is defined as 



(19) 



In the Majorana representation the charge counting term 
is then given by 

Tx{t) = iALi^[sin(A/2) a+r^f + cos(A/2) a+^f] , (20) 

where we introduced the new fields a± ~ {ai ± 02)/ V2- 
Because of the normalization of these new fields a fac- 
tor of \/2 has to be absorbed in every coupling constant. 



Analogously to the procedure shown in Ref. [431 we ex- 
press the adiabatic potential in terms of Green's functions 
(CPs) 



duj 
2^ 



dA 



D 



+ sin(A)L>„ - cos(A)i:>„ 



(21) 



where we used A — (A_ — A+)/2 and D represent ex- 
act GFs whereas g represent the GFs of zeroth order in 
coupling. Calculation of the exact GFs for the essentially 
quadratic Hamiltonian (we have set 7 = in the whole of 
this section) is straightforward but rather lengthy. The 
final result for the CGF reads 



00 



2i\ 



1) 



+ Ti(cj) [uRinp - 1) + npiuL - 1)] (e"'^ - l) 
+ Ti{uj) [nLinp - 1) + nrinn, - 1)] (e*^ - l) 

+ T2{uj)nR{nL - 1) (e"2*^ - l) }, (22) 

where Ti{uj) and T2{uj) are effective transmission coeffi- 
cients which are given by 



Ti{uj) = — and T2{uj) = — with (23) 
ao ao 



+ [KIj^ + KI]'^{kX[{Ai + A2f - 4cj2]2 + l6a;2(A2 + - 2ijj'^f) + 2bQ{A\ - Lo'^f{Al - Lo'^f, 
ai =2A-i^,{l6At(Ai + A2f{Lo^ - AiAa)^ + IQKWiAl + A^ - 2ujy + KXkI[{A^ + Aa)^ - Auj^Y], 
"2 =Klji{KX[{Ai + A2f - 4w2]2 + i6t^2j^2 + - 2lo'']^]. 

I 



In Eq. (|22|) the factors e**'*' and e^'^^^ distinguish between 
two different types of charge transfer processes: the first 
kind of processes transports a single charge across the 
system while the second one transports two charges. The 
"±" indicates the direction of charge transport where 
means "in direction of applied voltage" and "-" 
means "against direction of applied voltage". As ex- 
pected, one can see that in the zero temperature limit the 
terms corresponding to the "-" vanish as the Fermi func- 
tions become Heaviside's step functions. Unfortunately, 
at least in the chosen parameter range the FCS does not 
point towards collective processes involving three or four 
tunneling electrons. One reason for that might be the 



absence of mutual many-particle correlations between the 
dots4^i^ An inclusion of such interactions leads to the 
Toulouse point break down so that we have to postpone 
answering this question to a future work. 
We also note that in the case of Ai = = A2 our CGF, 
using the proper convention of notation, coincides with 
the result of [2g in the case of zero magnetic field. This 
means that without any magnetic field applied the double 
Kondo impurity behaves just as a single one. 
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A. Linear response regime V,T ^ 

In this case we can expand our FCS in the following 
way 

V 

lnx(A)«ry'^/HUo + ..., (24) 



where we write the integrand as a single function /(w) for 
simplicity. Thus the integration over uj becomes trivial 
if we limit ourselves to the zeroth expansion term. Then 
the FCS reads 

Inx(A) = ^In [1 + ri(c. = 0)(e'^ - 1) 

+ r2(c^ = 0)(e2^^-l)] , (25) 

which is equivalent to a manifestly binomial 
distribution^ 

X(A) = In [ 1 + ri(0)(e'^ - 1) + T2(0)(e2'^ - 1)] ^ 
= [1 + T,ie'^ - 1)] . (26) 

Te = yjTiiiiJ — 0) is an effective transmission coefficient. 
In terms of the introduced above the transmission co- 
efficients reduce to 

ao = [IGA^A^ + Kl{K\ + )(Ai + Aa)]^, (27) 
ai = 2(Ai + ^^fKlj^K\{\%^\t^\ 

+ iAi+A2fKlKl], (28) 

a2 = Ktj,Kl{A,+A2r. (29) 

At this point we would like to call the reader's atten- 
tion to the fact that for the case Ai = — A2 both Ti 
and T2 vanish whatever the choice of coupling constants. 
This means that in this case the transport through the 
system is suppressed and we observe an antiresonance at 
w = 0. This feature can be understood by the fact that 
for each electron which spin-flip tunnels across the sys- 
tem its corresponding holes does the same via the other 
quantum dot. This effect is very similar to the one found 
in the noninteracting double QDs.— (However, due to 
its different nature it does not allow to be used for spin 
filtering as suggested in Ref. [i^. See also Ref. [i^ for 
other interference effects.) Whether this feature pertains 
to the Toulouse point only is easily answered by a number 
of different perturbative expansions around this special 
point 1^ We explicitly performed the lowest order pertur- 
bative expansion in 7. This is a rather involved calcula- 
tion, details of which we present in the Appendix. In the 
case of small bias voltages the transmission is dominated 
by the lowest order in uj contributions from the respective 
self-energy. It turns out, that the only constant term is 
generated by the Toulouse point terms 7 = 0, higher or- 
der 7-terms coming with higher powers of uj. This can be 
understood as all Toulouse point correction terms being 
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FIG. 2: The nonlinear I{V) for the double QD setup at zero 
temperature. The plots show the current for Ai — 2 and 
A2 — 5. The different graphs correspond to K+ = and 
K- = 1 (dash-dotted line), K+ = 2 and K- = 2 (short 
dashed line), K+ = 1 and — (sohd line) and A'+ — 
and K- = (long dashed line). All energies are measured in 
units of isT^R. 

of inelastic origin kicking in at finite energies. This is, 
of course, consistent with the scaling dimensions of the 
correction terms. 

Very similar picture emerges upon loosening the re- 
strictions That is why we expect the antiresonance 
to be robust and universal beyond the Toulouse point. 
The only requirement for its realization is the fine-tuning 
of the magnetic fields to opposite values for both dots. 
We believe that this is experimentally feasible by e. g. 
applying inhomogeneous in-plane magnetic fields with fi- 
nite spacial gradient in the setups which were used in 
Refs. [HOill. 

B. Finite voltage, zero temperature 



Eq. ([22|) allows for the calculation of any cumulant de- 
sired for the given setup. At this point we would like to 
present some experimentally observable quantities and 
discuss their specific features. In Fig. [2] the electric cur- 
rent through the system is shown for four different sets 
of coupling constants while the strength of the magnetic 
field is the same for all plots. The basic behaviour of all 
four curves is qualitatively identical: the current shows 
up two distinct plateaus as well as two steep shoulders 
before it saturates for a sufficiently large voltage. In 
any case the two voltages where the current increases are 
equal to the strength of either of the magnetic fields. The 
reason for that kind of behavior is rather simple. Since 
the primary transport mechanism is the spin-flip tunnel- 
ing it is suppressed for both dots at voltages smaller than 
the respective magnetic field. As soon as V overcomes 
the smaller of the fields, the respective transport channel 
opens and there is a rapid increase of the current. The 
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F(V) F(V) 




FIG. 3: Fano factor for fixed magnetic fields Ai — 2 and 
A2 = 5 as a function of bias voltage. The different graphs 
correspond to K+ = = if- (dash-dotted line), K+ = 0.5 = 
K- (short dashed line), K+ = 0.5 and K- = (solid line) 
and K+ = and K- — 0.5 (long dashed line). Solid line 
F{V) = 0.5 is guide to the eye only. 



second step is then associated with the other field. The 
overall saturation of the current is related to the finite 
total spectral density of the constriction as required by 
the sum rules. 

As far as the parameter dependence is concerned, see 
Fig. [21 the general trend that can be observed is that the 
shoulders of the curves tend to smear out for larger val- 
ues of if+ and K- (all couplings are measured in units 
of Klu, parameters with the dimension of an energy are 
measured in units of Kj^j^). Both step- like features are 
most pronounced in the case of K± = 0, which is the 
case of only spin-flip tunneling term present in the orig- 
inal Hamiltonian. A finite K± corresponds to additional 
transversal coupling of the impurity spins, which induces 
spin precession in the free case and is a source of inde- 
pendent spin-flips in the coupled case. This effect washes 
out the steps in the I — V characteristics. Interestingly, 
only one of K± being finite seems to cause a significantly 
weaker "dcphasing" than a situation of both of them be- 
ing non-zero. 



FIG. 4: Main graph: Fano factor as a function of bias 
voltage for different temperatures and fixed magnetic fields 
Ai = 2 and A2 = 5, T = 0.005 (dash-dotted line), T = 0.05 
(solid line), T = O.f (dashed line). The divergence of the 
Fano factor for V can be explained by expanding noise 
and current (see text). Inset: Current as a function of bias 
voltage at different temperatures: T = 0.005 (dash-dotted 
line), T = 0.5 (solid line) and T — 1. (dashed line). Solid line 
F{V) = 0.5 is guide to the eye only. 



observed in the transport through constrictions with in- 
ternal degrees of freedom, e. g. it is known to appear in 
the Fano factor of the resonant level setupi^ 



C. Finite temperature effects 

To access the Johnson-Nyquist noise we set V = and 
assume T to be small.— 1^ In this case the FCS reads 



lnx(A)= / — In l + nFil-np) 



x[T2(w)(e 





2iX 



-2i\ 



-2) + 2Ti{u){e'^ + e 



Now we use the fact that np{l — np) = —p^^d^np {(3 = 
1 /T is the inverse temperature) and calculate the noise. 
Partial integration of the resulting expression gives 



The Fano factor F{V) as defined in Eq. (US]) and I{V) 
show up similar features at V = A 1.2, sec Fig. [31 

For small voltages and only spin-flip processes present 
K± = the Fano factor approaches unity. For finite 
transversal coupling K± ^ this value becomes nonuni- 
versal. Interestingly, the effect of finite Kj^ is much more 
pronounced than that of The reason for that is 

the fact that while finite Kj^ indicates the presence of 
the transversal couplings to the individual terminals, A'-_ 
measures its asymmetry. 

On the contrary, at I^ — )• oo the Fano factor reaches 
the universal asymptotic value 1/2, whatever the cou- 
pling strengths and temperature. This effect is usually 



S « 4/3-i-[Ti(0) -I- 2r2(0)] = , (30) 

where we identify ri(0) + 2r2(0) with the effective trans- 
mission coefficient from Section UlIAl 

In Fig. [H we show the behaviour of the Fano factor as 
a function of voltage for different temperatures. As ex- 
pected the Fano factor grows for increasing temperatures 
due to the onset of thermal fluctuations. The tempera- 
tures chosen are much smaller than if la, which repre- 
sents the Kondo temperature because any T comparable 
or even higher than if^i? eliminates all features from the 
plot as already can be seen from the inset. In addition to 
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the features aheady discussed in the T = case the Fano 
factor rises for V 0. This behaviour can be explained 
by the fact that the expansion of the current in V starts 
with a second order term while the noise starts with a 
first order term. 



IV. CONCLUSIONS 

To conclude, we have discussed the nonlinear trans- 
port properties of a double quantum dot in the Kondo 
regime. Using a series of dedicated transformations we 
rewrote the original Hamiltonian in one special region of 
the parameter space in terms of a quadratic Majorana 
resonant level model, which is conveniently diagonaliz- 
able even under nonequilibrium conditions. We have 
explicitly calculated the generating function of the full 
counting statistics and discussed its most prominent fea- 
tures, which reveal themselves in individual cumulants of 
the charge transport. In particular, we find a full sup- 
pression of transport for a special constellation of applied 
magnetic fields. We argue that this antiresonance feature 
is robust and universal even beyond the Toulouse point 
by an explicit perturbative expansion around it. One 
possible route for further progress would be an analy- 
sis of this phenomenon by numerical and possibly more 
advanced analytical methods for the generic system pa- 
rameters (e. g. by quantum Monte Carlo methods or 
functional renormalization group approaches) . 
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Appendix 

Here we want to summarize the most important details 
of the perturbative expansion around the Toulouse point. 
The full fledged calculation of the corresponding I{V) 
turns out to be very complex so that we employ the fol- 
lowing approximation strategy. From the Toulouse point 
calculation we know that (at least at zero temperature) 
the electric current through the system is given by an 
energy integral over the voltage window of the imaginary 
part of the QD retarded GF 



most conveniently done in the Keldysh formalism from 
the following GFs, 



i{Tc a+(t)a+(t') exp | - 7 dh 



D—'-+ - 

a-i-Q-i- 

[:^l(ti)Vs(<i) : [a+(ti)6+(ii) + a_(ti)6_(ii)]]}), 

(32) 

which we expand for small 7. The first order vanishes 
due to the normal ordering in the free fermion sector. 
The second order we find to be given by 

■ 2 /• 

{Tca+{t)a+{t') [a+16+1 + a_i6_i] [a+2&+2 + 0-2^-2]), 

where the second index is the label of the time argu- 
ment. Schematically it corresponds to a simple com- 
pound Majorana-free fermion loop diagrami^ Next we 
use Wick's theorem and find the first expectation value 
of the free fermion operators to be the product of two 
GF. The Majorana part is a bit more involved. Since 
many contributions vanish for Ai = — A2 one finds a 
manageable expression 

- iDa+a+ {t - h)[Db_a^ {1-2 - t')Db^a^ {h - (33) 

-Db^b^{h^t2)Da^a+{t2-t')] (34) 

+iDa+b^ {t - ti)[Db_a+ {t2 - t')Da_a^ {h - t2) (35) 

-Da+a+{t2~t')Da_b^(ti~t2)\, (36) 

where all GFs are to be understood as matrices in the 
Keldysh space. Transformed into Fourier space we are 
left with expressions of the type 

Diiuj) j deD2{e) j dnG^{uj-€ + n)G^{n)D^{uj) , (37) 

where are the (local, taken at x = 0) free fermion GFs 
of i/^s-fields, which are known, see e. g. Rcf. [l^. D1.2 
are the "outer" GFs from Eq. ([36|) . which correspond to 
the one with the time arguments t — ti and t2 — t' in the 
time domain, while D2 represents the "inner" GF with 
the argument ti — ^2 • 

The calculation of the dot GFs is somewhat lengthy 
but straightforward by writing down the action 



5 = %,?] + /^[E"K".'(^+'-) 



-iA_{a^h+ + a+b^)+iK+i-jssb++iKLR£.fa+ , (38) 



^?+a+ = -*')({«+ W> «+(*')}) , (31) where we defined A± = Ai ± A2 and the superfield 

(ai_, 0;+, 6i+), and integrating out the lead 



which plays the role of the effective transmission coeffi- 
cient (of course, this is compatible to the Wingreen-Meir 
formula^). So we just need corrections to this GF. It is 



fermions. ^^(((A+jw) denotes the 4 x 4-matrix GF for 
the constriction without couplings to the electrodes. In- 
tegrating them out we can read off the necessary GFs 
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from the GF for the superfields which has the structure: 



D 



^ Ddj^dj^ Da_\-h^ I-^a^a- ^aj^b- \ 

Db+a+ Di^i^ Di^a- 

Da._ Q-i- 6-1- Da_b- 

\Dh_aj, Dh-a- ^b-b-J 



(39) 



where we indicated those components that are zero with a 
hat. Next we use the above GF and expand the structures 
(|57|) to the fowest order in w. In total one needs to take 

care of 16 different terms of the type ((37|) for both 

and H — V component at w — > 0. AU but two terms yield 
higher order contributions to the expansion. The only 
term in the component is the expression 



-{n-e)G—{n). (40) 



Performing the 51-integration first leaves us with an in- 
tegral over -D(7Ja_(^)- This GF is an odd function in e 



and has to be multiplied with a function that includes 
the bandwidth but is even in e. This means that the in- 
tegral vanishes since integration and taking the limit for 
the bandwidth can be exchanged. 

The only term that is non-zero in the lowest order in 
Lo is given by 



^D-;,^D+Xj deD-+^_{e) J dnG-+ - e)G+- (n) 



de £3 



(A2 



4 I ,oA/5^ 



const + O(y^). (41) 



That means that this contribution vanishes for small volt- 
ages much faster than the leading order term and there- 
fore ensures the existence of the perfect antircsonance in 
vicinity of the Toulouse point. 
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